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ABSTRACT 

We present the results from a multi-epoch survey of two regions of M33 using the 3.5m 
WIYN telescope. The inner field is located close to the centre of the galaxy, with the 
outer region situated about 5.1 kpc away in the southern spiral arm, allowing us to 
sample a large metallicity range. We have data for 167 fundamental mode Cepheids 
in the two regions. The reddening-free Wesenheit magnitude Wm period-luminosity 
relations were used to establish the distance modulus of each region, with pinner = 
24.37 ± 0.02 mag and pouter = 24.54 ± 0.03 mag. The apparent discrepancy between 
these two results can be explained by the significant metallicity gradient of the galaxy. 
We determine a value for the metallicity parameter of the Period-Luminosity relation 

7 = ^ — = —0.29 ± 0.11 mag dex~^, consistent with previous measurements. 

dlogZ 

This leads to a metallicity corrected distance modulus to M33 of /i^ = 24.53 ± 0.11 
mag. 

Key words: catalogues ~ Cepheids - galaxies: individual: M33 - galaxies: distances 
and redshifts 



1 INTRODUCTION 

Cepheid variables have been used to determine gal actic dis- 
tances for many decades (e.g. lFerni^ |l969), Frccdm an et al.l 
l|200lll and references therein). The relationship between 
their period of pulsation and luminosity has been consid- 
ered to be well defined at several times during its lifetime. 
However, the effect of the Cepheid's chemical composition 
on the Period-Luminosity (PL) relation has never been ade- 
quately determined, leaving the relation open to systematic 
uncertainties that have not been fully investigated. Proper 
calibration of this relation is essential not just to explain 
how metallicity affects the pulsations of Cepheids, but to tie 
down the zero point of the extra-galactic distance scale and 
thus properly measure the Hubble constant. 

Metallicity and differential absorption will both change 
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stell ar colours |Bell fc Rodeerj Il969l ) . Theoretical models 
(e.g. ICaputo et al.l l|2000l )) show that the metal content of a 
Cepheid will affect the derived distance modulus. If a star 
of an unknown metallicity is at an undetermined distance 
then it is difficult to untangle these two effects on the stars 
colours. Ideally the PL relation would remove this problem 
as the distance to the star could be found without calcu- 
lating the colour of the Cepheid. The problem arises as the 
calibration of the PL relation does not take this effect into 
account. The distances that were used to calibrate the zero 
point of the relation were determined without considering 
how the reddening of the Cepheids was affected by the com- 
bination of these two mechanisms, meaning that a system- 
atic offset may have been introduced into the relation. 

There has been some investigation into the effect 
that metallicity has o n the PL relation, be gi nning with 
Freedman fc Madord (|l990l ). iBeaulieu etahl (|l997l ) and 
Sasselov et al.l ( 19971 ) used the Cepheids detected in the 
EROS gravitational lensing survey to show that discrepant 
values of Hq could be reconciled by appl ying a metallic- 
ity co rrection to Cepheid distance moduli. Kennicutt et al.l 
(|l998h used the Hubble Space Telescope to observe Cepheids 
in M 101, a nearby spiral galaxy with active star forma- 
tion regions but low dust extinction properties. Choosing 
a galaxy with low extinction means that there would be 
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less reddening from dust; any reddening of the colour of the 
Cepheids could be explained as a metallicity effect. This 
study found that 7 — —0.24 ± 0.16 mag dex~^ (see also 
Freedman et al., 2001). However, their conclusions were 
based on small n umbers of Cepheids , combined with the 
PL relations from lUdalski et al.l (|l999l ). It is far from ideal 
to use such a small sample to quantify an effect that could 
mean that the entire extra-galactic distance scale is sub- 
ject to a systematic offset. Usi ng a similar meth od but 
with a larger sample of Cepheids, iMacri et al.l l|2006l ') found 
7 — —0.29 ±0.09 mag dex~^. Using a la r ge sam ple of nearly 
700 Cepheids in 17 galaxies, iKochanekl (119971 ) had found a 
result compatible with these other works. 

Most works use a sample of Cepheids covering a range in 
metallicity and e nforce consistency between di stance mod- 
uli; th is is what iKennicutt et al.[ l|l998l ') and iMacri et al.l 
ll2006t) did for instan ce. An alternative approach used e.g. 
by Sakai et al.l (|2004l ) is to assume that the distances to sev- 
eral galaxies are known a priori, from the tip of the red giant 
branch in their case, and find the correction needed to bring 
Cepheid distances to these galaxies in agreement with that 
other standard candle. 

There are two possible solutions to the problem of 
the metallicity-extinction d egeneracy. The method followed 
by IKennicutt et all (|1998D is one possibility, the elimina- 
tion of extinction from the equation. Theoretically this 
method should work, but in practice it is not so simple. 
IKennicutt et al.l (|l998l ) chose a galaxy with low extinction, 
not zero extinction. This implies that the extinction is not 
entirely eliminated, leaving the degeneracy that we are try- 
ing to overcome. However, this degeneracy could be removed 
by using the reddening-free Wesenheit magnitude Wvi PL 
relation to calculate the distance. By eliminating extinction 
in this way, any discrepancies in the apparent distance mod- 
ulus could be explained by the metallicity effect. 

The work presented here is based on a multi-epoch sur- 
vey of M33 conducted with the WIYN 3.5 m telescope. M33 
is a well studied galaxy, with observ ations of Ceph eids per- 
formed by Hubble as early as 1926 (lHubblelll926l '). As it is 
nearby, many studies have determined the distance of M33 
using differe nt and independent m ethods, such as eclips- 
ing binaries (|Bona nos et al. 20061). the tip of the red gi- 
ant branch (TRGB ) (Broo ks et al.ll2004 ICalleti et al.ll2004 
iMcConachie et all 120041 ). masers JBrunthaler et all |2005| ) 
and Cepheids themselves (|Freedman et al.lll99ll ). Although 
some of these methods may not provide the same accuracy 
that a Cepheid distance does, they allow us to have an inde- 
pendent starting point for our analysis of the PL relation of 
M33. One of the reasons that M33 was chosen for this study 
is that it is well studied. Its inclination angle is 53° , and so 
any effect on distance modulus (/i) from the orientation will 
be extremely small. The effect on /i due to the orientation of 
the galaxy was found to be very small (less than 0.01 mag). 
The metallicity gradient across the galaxy is well defined, so 
the chemical composition of the Cepheids is relatively easy 
to determine just from their position in the galaxy. These 
properties inherent to M33 make it the ideal target for a 
study of this nature. 

The paper is organised as follows. Section [2] describes 
our observations. Section |3] describes the reduction process, 
including the development of the photometry pipeline and 
the final calibration and astrometry. Section |4] contains the 




Figure 1. The two regions observed for this survey^. The inner 
region contains fields 1 and 2; the outer region contains fields 
3 and 4. The fields were chosen such that they would cover the 
whole metallicity range of the galaxy. Orientation: North up. East 
right. 



final catalogues and assesses the accuracy of the photometry, 
with the final PL relations and distances in section[5l Finally, 
we examine the effect of metallicity on our results and derive 
a value for the metallicity correction parameter 7. 



2 OBSERVATIONS 

The data set consists of 25 epochs of BVIc images of four 
fields of M33 obtained between 1998 and 2001 on the 3.5 
m WIYN telescope. The first 20 epochs used the WIYN 
s2kb Imager, a single thinned Tek/STIS 2048x2048 CCD 
with 21 /im pixels. For the final 5 epochs the new Mini- 
Mosaic (MiMo) camera was used. The Mini-Mosaic consists 
of two thinned SITe 2048x4096 CCDs with 15 /im pixels. A 
single exposure covered 6.8 x 6.8 arcmin for the Imager, and 
9.6 X 9.6 arcmin for the Mini-Mosaic, with pixel scales of 
0.195 and 0.141 arcsec pixel"^ respectively. 

Four fields were observed for the survey, two in the cen- 
tre of M33, and two in the outer, southern part of the galaxy. 
The locations of the fields are shown in Figure[l] The J2000.0 
coordinates of the fields are shown in Table [1] All the fields 
were observed in the B, V and Ic filters with both cam- 
eras. The position of the fields were chosen such that the 
metal-rich population at the centre of the galax;y could be 
compared with the more metal-poor population in the outer 
region. Cepheids tend to be concentrated in the spiral arms, 
so the outer fields were located in the southern part of the 
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Table 1. J2000.0 coordinates of the four fields. M33-1 and M33-2 
are situated in the centre of the galaxy. M33-3 and M33-4 are in 
the southern spiral arm. Deprojected angular and linear distances 
of each field to the centre of the galaxy are given, assuming M33 
is at a distance of 840 kpc, has an inclination angle of 53° and 
the major axis has a position angle of 22°. 



Field 


a 


5 


Angular 


Linear 


Imager 


M33-1 
M33-2 
M33-3 
M33-4 


l'"33'"49.12* 
l'"34™06.13= 
l'^33™48.99= 
l''33'"25.63'' 


30° 42' 56. 48" 
30°38'52.50" 
30°20'24.58" 
30°20'25.25" 


217.64" 
217.10" 
1242.13" 
1290.30" 


0.886 kpc 
0.884 kpc 
5.059 kpc 
5.255 kpc 


Mi-Mo 


M33-1 
M33-2 
M33-3 
M33-4 


l'"33'"51.77^ 
1*^ 34"" 16.75* 
l'*34'"00.16= 
l'*33™21.13= 


30° 45' 16.56" 
30°36'35.64" 
30°21'46.08" 
30°21'52.20" 


367.50" 
409.13" 
1161.16" 
1220.28" 


1.497 kpc 
1.666 kpc 
4.729 kpc 
4.970 kpc 



Table 2. Log of Observations. Table [2] is available in full elec- 
tronically. Cameras: IM = S2KB Imager, MiMo = Mini-Mosaic 



UT Date 




MHJD 


Field 


Filter 


Cam. 


Seeing 


1998 Aug 


19 


51044.427 


M33-1 


B 


IM 


0.53" 


1998 Aug 


19 


51044.433 


M33-1 


Ic 


IM 


0.53" 



galaxy to coincide with the clearest arm. The positions of 
the fields and the addition of the Mini- Mosaic camera meant 
that some objects were observed in either both inner or both 
outer fields, giving up to 38 observations for objects in the 
overlap regions. 



3 DATA REDUCTION 

Preliminary reduction of the data were performed using 
the routines in the lRAl{3 CCDPROC package. This con- 
sisted of overscan correction, bias subtraction and flat field- 
ing. The p hotomet r y wa s performed using the standalone 
DAOPHOT (jStetsonl Il987l ) and allframe (jStetsonl Il994l ) 
packages kindly provided by P. Stetson. 



3.1 Photometry Pipeline 

The data set used in this project is fairly large, containing 
over 300 images of the four fields. To process such a large 
data set a completely automatic photometry pipeline was 
developed, covering all the steps from initial star selection 
to the final PSF photometry. 

Initial star lists were compiled using the find routine in 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



DAOPHOT with a 3(T detection limit. Aperture photometry 
was performed on each list to get rough magnitudes using 
an 11 pixel aperture and sky values from the annulus from 
19 to 23 pixels. 

The creation of the point spread function (PSF) model 
is the most important part of the process. The fields are 
extremely crowded; the inner field has an average density of 
890 stars arcmin"^, or approximately 10 pixels between each 
centroid position. Of course these calculations do not take 
into account the dumpiness of the distribution, the galactic 
centre has a much higher stellar density than further out. 
As the fields are so crowded it is imperative that the PSF 
model is as accurate a representation of the stellar profile as 
possible, but it also means that it becomes a more taxing 
problem. The stars chosen to make the model should be 
the brightest, least crowded available, and with such heavy 
crowding these are very rare objects. 

The problem was overcome using a combination of the 
DAOPHOT PICK and PSF routines along with a series of scripts 
to analyse the resulting fits. The initial model was created 
by allowing pick to choose the 200 best stars from the aper- 
ture photometry file then using PSF to make six model pro- 
files using the difi^erent Gaussian, Moffat and Lorentz func- 
tions available, allowing the model to vary linearly across the 
frames. The residuals of the fit were analysed and any stars 
fiagged as saturated or with fitting residuals more than twice 
the root-mean-squared (RMS) residual value were removed. 
All the stars except the PSF stars were then removed from 
the images and the PSF was remade with the new star list 
from the cleaned image. The process of analysing the resid- 
uals, removing neighbour stars and remaking the PSF was 
repeated a further five times and resulted in a clean, well 
fitting model. Allowing the model to vary quadratically and 
repeating the cleaning process more times were both tested 
but did not significantly improve the results. 

Profile fitting photometry was performed on each frame 
separately using the ALLSTAR routine, which produced an 
image with all the fitted stellar profiles subtracted. A fitting 
radius of 0.78 arcsec was used, corresponding to the typical 
seeing. This new image was run through the star finding pro- 
cess again to detect any stars that were too faint or crowded 
to be seen in the first pass. This new star list was appended 
to the original and run through ALLSTAR again. 

Once reasonable PSF magnitudes had been obtained 
for each frame, the images were combined using montage2 
(Stetson, personal communication) using the best 50 frames 
to produce two master images, one for the inner and another 
for the outer fields. These deep, medianed images were then 
run through the pipeline described above to produce mas- 
ter star lists for the two regions. The master lists contained 
over 200,000 stars in the inner region and around 100,000 for 
the outer region. The lists produced from these frames were 
used rather than the original lists for two reasons. Firstly, a 
medianed image made from many frames has a much better 
signal to noise ratio (S/N) than a single frame. Secondly, us- 
ing identical star lists for each frame and each filter is nec- 
essary when using the allframe package, allframe was 
used to calculate the final PSF magnitudes from 'fixed po- 
sition' coordinate lists. 

The final step of the instrumental photometry was to 
calculate the aperture corrections and set all the photome- 
try to the same zero point. Landolt standards were observed 
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on 1999 October 3rd so this night was chosen to be the ref- 
erence night. Aperture corrections were calculated for each 
frame take n on the refere nce night (4 fields in 3 filters) using 
DAOGROW (|Stetsonlll99Cll '). Using the ALLFRAME photometry 
files as input, all the stars in each frame except the stars used 
to make the PSF model were subtracted. The PSF stars are 
ideal for calculating the aperture correction as they are by 
their nature bright and relatively isolated. The stars remain- 
ing in the frame were then measured in nine apertures with 
radii between 4 and 20 pixels, with the sky measured from an 
annulus between 25 and 30 p ixels. Th e output photometry 
was then put into DAOGROW (jStetson 1990) to calculate the 
growth curves and derive the total magnitudes. The differ- 
ence between the total magnitude and the PSF magnitude 
for each star was calculated, and any bad outliers were re- 
moved. These stars could have been affected by nearby stars 
that had been missed, or even nearby PSF stars in the large 
aperture. The best value of the correction was found using a 
least-squares fit. The corrections were then subtracted from 
the PSF magnitude to put the ALLFRAME photometry on 
the same scale as the aperture magnitudes. The aperture 
corrections for the reference night had an average value of 
0.035 mag. 

To get correct photometry for every night the zero-point 
offset between each frame and its corresponding reference 
frame was calculated. As the fields overlap, field 1 was cho- 
sen as the reference for the inner region and field 3 for the 
outer region. This ensured that stars that were observed in 
different fields on the same night were on the same zero 
point. The zero-point offset of each frame was found by tak- 
ing the average difference between each magnitude on the 
reference frame and the frame in question, typically using 
1000-2000 stars. The star furthest away from the average 
was removed and the average recalculated until the disper- 
sion of the differences was less than 0.01 mag. This was 
found to be a very robust method, distinguishing zero-point 
differences down to one hundredth of a magnitude. 

3.2 Photometric Calibration and Astrometry 

Observations of standard fields l|LandoltllT993 ) were taken 
on 1999 October 3rd with the Imager CCD. Three standard 
fields were observed repeatedly, with a total of 17 different 
standard stars. Aperture photometry was performed on the 
standards with the same settings as the aperture correction, 
that is a 20 pixel aperture and a 25 to 30 pixel sky annu- 
lus. The PHOTCAL routines in iraf were used to construct 
standard solutions of the form; 



Table 3. Photometric solution of 1999 October 3rd for the Imager 
CCD 



rricai + a + b X Xm + cx Colour + d x UT (1) 



where niinst is the instrumental magnitude, rricai is the cal- 
ibrated magnitude, Xm is the airmass and Colour corre- 
sponds to [B — V) for the B and V frames and (V — Ic) 
for Ic ■ The coefficients for each transformation equation are 
shown in Table [S] The solution was applied to all Imager 
frames using invertfit. The UT term was included in the 
solution as the extinction changed very slowly over the ref- 
erence night. Although the coefficient of this term in the 
solution is much smaller than the others we believe that it 
is still important to include it. The aim of the process is 
to get the most accurate photometry possible, so if an ef- 
fect is present and can be measured on a scale comparable 



Filter 


a 


b 


c 


d 


RMS 


B 


24.1552 


0.2760 


-0.0682 


-0.0062 


0.0098 


V 


24.3051 


0.2006 


-0.0123 


-0.0062 


0.0107 


Ic 


23.7141 


0.0903 


-0.0545 


-0.0026 


0.0187 



Table 4. Photometric solution of 2000 October 3rd for the Mini- 
Mosaic Camera 



Filter 


a 


b 


c 


RMS 


Chip 1 


B 


23.9118 


-0.0935 


0.0172 


0.0494 


V 


24.1484 


-0.0235 


0.0079 


0.0359 


Ic 


23.6443 


-0.0374 


-0.0083 


0.0472 


Chip 2 


B 


23.9120 


-0.0587 


-0.0143 


0.0486 


V 


24.1465 


0.0128 


-0.0245 


0.0302 


Ic 


23.6373 


-0.0610 


0.0042 


0.0455 



to the accuracy we are trying to achieve then it should be 
considered. 

This calibration is only appropriate for the frames taken 
with the Imager CCD. Once applied to the Imager frames, 
the stars on the reference frames could be used as local stan- 
dards to derive a photometric solution for the two Mini- 
Mosaic CCDs. The reference night for the Mini-Mosaic was 
chosen to be 2000 October 3rd as this had the best image 
quality for all fields. Stars that appeared in all frames of field 
3 on both the photometric reference night and on the Mo- 
saic reference night with ay < 0.02 mag were selected as the 
local standards - around 200 stars for each chip. The PHOT- 
CAL routines were used to find the solution of the transfor- 
mation equations from the Mosaic instrumental magnitudes 
to the Imager Landolt calibration. For this transformation 
the time dependence and airmass term are absorbed into the 
zero-point coefficient of this solution and a quadratic colour 
term was added, resulting in transformation equations of the 
form 

rriMiMo ~ mim + a + b x Colour + c x (Colour)'^ (2) 

where rriMiMo refers to the magnitude on the Mini-Mosaic 
instrumental scale and mim refers to the calibrated Imager 
scale. The quadratic colour term of the Mini - Mosaic camer a 
is a feature that has been described before by 'stetson' (|2005h . 
The coefficients of the equations are shown in Table |3] 

The stacked frames created for making the ALLFRAME 
master lists were also used to perform the astrometry. For 
the outer field the USNO-B catalogue (|Monet et al.ll2003l 'l 
were used, resulting in a solution with RMS values of 0.4". 
The USNO-B catalogue was not appropriate for the inner 
field as it covers the central part of the galaxy, which is 
extremely crowded. For this region the astrometric solution 
w as found by m a tching by eye to the catalogues published 
in iMassev et all l|2006l ). They also use USNO-B catalogue 
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Table 6 . Typical magnitudes and uncertainties from the weighted 
averages catalogue. 



Field 


Filter 


a = 1% 


(T = 10% 


Inner 


B 


21.125 


24.129 




V 


21.693 


23.854 




Ic 


19.896 


22.597 


Outer 


B 


21.336 


24.404 




V 


22.136 


24.195 




Ic 


20.840 


22.643 



for this region but have many fewer stars. This process gave 
a much better solution, with an RMS below 0.07". 



4 RESULTS 

4.1 The Catalogues 

Two catalogues are available from the WIYN M33 Survey. 
The first contains weighted average B, V and Ic magnitudes 
for each star. A sample of this catalogue is shown in Tablet^. 
The second contains all the individual measurements for all 
the variables cross-id entified with the Dee p CFHT Photo- 
metric Survey of M33 l|Hartman et al.|[2006l ). The catalogues 
are available in full electronically. They cover the magnitude 
range between 13.8 <V< 26.7, with 75465 and 146820 in 
the final outer and inner field lists respectively. The magni- 
tude range of the inner field is limited slightly by the degree 
of crowding, and only contains objects down to V = 25.7. 
Colour-magnitude diagrams for both fields are shown in Fig- 
ure [1 

The uncertainties given in the average photometry cat- 
alogue are the errors on the weighted average magnitudes, 
calculated from the allframe results. Figure |3] shows the 
relation between magnitude and uncertainty for the V band 
in the outer field. A second series becomes visible around 
mv = 20 mag. These objects have fewer measurements, a 
maximum of five, as they appear only on the Mini-Mosaic 
frames (which has a large field of view than the Imager). 
However, a comparison between magnitudes calculated from 
the Imager and the Mini-Mosaic shows the photometry is 
compatible. Neglecting the second series we obtain 1% pho- 
tometry for objects brighter than my ~ 22, and 10% down 
to mv ~ 24. The full analysis is shown in Table [S] These 
results are compatible with theoretical estimates from signal 
to noise calculations. 

4.2 Comparison with others 

To check the accuracy of our results the a verage magni- 
tudes were compared with those obtained by iMassev et al.l 
l|2006t ). Their survey covers a larger area than the survey 
presented here but does not go as deep. The comparison 
for each filter in both fields is shown in Figure 3] The 
photometry from this work is considered as the reference, 
Amag = mscomcroft — rnMasaey These plots show that the 
photometry is sound. Most objects have a difference near 

^ The full catalogues are available on request from the authors. 




15 EO S5 



Figure 3. Magnitude - a relation for V in the outer field. The 
second series with slightly higher uncertainties for a given mag- 
nitude is due to the Mini-Mosaic magnitudes dominating these 
objects, as these stars have few measurements. 



Table 7. Differen ces between weighte d average magnitudes in 
this catalogue and iMassev et al.l | |2006|) . Amag was found from 
a least squares fit to the brightest end of the comparison, the 
range is given in column 4. CAmag is the uncertainty on Amag. 
Outliers with Amag > ±0.5 were rejected, as they are most likely 
misidentifications due to the crowding in the inner field. 



Field 


Filter 


Amag 




Range 


Inner 


B 


0.008 


0.005 


16-19 




V 


-0.009 


0.008 


15-18 




Ic 


0.031 


0.008 


15-18 


Outer 


B 


-0.013 


0.008 


16-19 




V 


-0.014 


0.005 


16-19 




Ic 


0.021 


0.005 


16-19 



zero as expected. For the inner field in particular, there 
seems to be a secondary sequence near Amag ^ —0.8. We 
attribute this to misidentifications due to the high level of 
crowding. By excluding the central 4 arcmin of the galaxy 
from the comparison the effect was significantly reduced. 

The average differences were found using a least squares 
fit to the brightest three magnitudes in the comparison and 
are shown in Table [T] Only the brightest stars in the com- 
parison were chosen as the two catalogues are optimised 
for different magnitude ranges; were we to consider the 
whole r ange we would be b asing the fit on the very faintest 
stars in lMassev et al.l (|2006| ). However, Amag gradients were 
checked for, and the differences for each field-filter combi- 
nation were found to be constant over the whole range of 
magnitudes. We find a larger average Amag for Ic which 
is believed to be because the Ic frames are more crowded 
than B and V , particularly in the inner field. 
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Table 5. Sample of the weighted averages catalogue. The uncertainties are from the weighted average calculations derived 
from the ALLFRAME results. Where no observations were made the magnitude and error are represented by 99.999 and 9.999 
respectively. 



ID 


a hh:mm:ss 


5 dd:mm:ss 


B 




V 


(7V 


Ic 




Nb 


Nv 


Ni 


OUT-14501 


1:33:10.21 


30:17:22.7 


24.946 


0.423 


23.733 


0.161 


99.999 


9.999 


1 


1 





OUT-14510 


1:33:10.25 


30:17:05.8 


22.484 


0.155 


22.269 


0.096 


99.999 


9.999 


2 


2 





OUT-14523 


1:33:10.26 


30:17:07.2 


22.525 


0.092 


22.872 


0.075 


99.999 


9.999 


2 


2 





OUT-14527 


1:33:10.19 


30:17:50.4 


23.608 


0.312 


24.306 


0.248 


99.999 


9.999 


1 


1 





OUT-14532 


1:33:10.26 


30:17:15.2 


23.475 


0.109 


22.943 


0.067 


99.999 


9.999 


2 


2 





OUT-14539 


1:33:10.27 


30:17:17.2 


23.474 


0.125 


23.576 


0.098 


23.232 


0.335 


2 


2 


1 


OUT-14546 
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Figure 2. Hess diagrams for the inner (left) and outer (right) fields, with 0.025 mag bins. The inner field does not go as deep as the 
outer field as the degree of crowding is much higher. The effects of crowding and differential reddening are obvious in the inner field but 
are much less so in the outer field 



5 CEPHEID PERIOD-LUMINOSITY 
RELATION 

The aim of the survey was to detect Cepheids in the two 
fields and assess the effect that metalhcity has on the PL 
relation. The positions of the two fields were chosen such 
that the inner field covered the metal-rich centre, which has 
an [Fe/H] value slightly larger than the Milky Way (MW), 
and the outer field covered the southern spiral arm where 
the metallicity drops to slightly below that of the Large 
Magellanic Cloud (LMC). 



5.1 Cepheid Selection 

Cepheids were identified by cross referencing the catalogue 
with the fundamental mode Cepheid cata logue produced 
as pa rt of the deep CFHT survey of M33 (|Hartman et al.l 
I2OO6I ). The fundamental mode Cepheids we re identified us- 
ing Fourier parameters (jSimon First overtone 



Cepheids were excluded from the analysis as they are known 
to sit slightly higher in the log P — Mag plane than funda- 
mental Cepheids and would bias any zero-point fitting to- 
wards brighter magnitudes. 

To ensure an accurate PL relation the average magni- 
tudes for each Cepheid were recalculated by fitting Fourier 
series to the light curves, using the form shown in Equation 
[S] where rn{(j)) is the magnitude at phase (j) and a becomes 
the phase averaged magnitude. 

3 

m{(j)) — a + '^^hiCos{2-Ki(j) + Ci) (3) 

The light curve of each Cepheid was examined visually to 
remove any misidentifications or ones with incorrect periods. 
The amplitude of each Cepheid was calculated and Cepheids 
with Ab > Av > Ai^ were incl uded in the final s ample , 
similar to the selection criteria of iPietrzvnski et all (I2OO2D . 
Amplitude ratio selection is an important step as large num- 
bers of classical Cepheids are known to exist in binary sys- 
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Figure 4. Comparison between this work and the catalogues of iMassev et alj CO OG') for the inner field (left) and the outer field (right). 
This worlc is considered to be the reference and is denoted by Mg . iMasseT'etaL , (2006) by Mm and Amag = Mg — Mm- The average 
number of stars compared is 9700 for the inner field and 5100 for the outer field, corresponding to roughly 5% of our catalogues. All 
matched stars with a < 0.5 mag that were outside the central 4 arcmin were included. The secondary Ms — Mm < trends that are visible 
in the left hand panel are believed to be due to the severe crowding of the region, as the matching algorithm used will tend to match to 
slightly brighter stars in highly crowded regions. 



tems (|Szabadosll200l '). Very close companions would not be 
identified as separate stars in the initial photometry. They 
can however be identified by their characteristic effect on the 
Cepheid light curve. A Cepheid with an unidentified com- 
panion will still appear variable, but the hght curve would 
not reach the minimum value expected. The effect will be 
seen to different extents in the different photometric bands 
depending on the colour of the companion. This explanation 
can also be extended to blending; if a Cepheid is within a 
background of unidentified faint stars these objects will af- 
fect the shape of the hght curve in the same way. By requir- 
ing that the amplitudes of the Cepheids follow the criteria 
set out above we can remove these effects from the sample. 
The visual inspection and amplitude selection left 91 and 28 
Cepheids in the inner and outer field samples respectively. 



5.2 Metallicity effects 

Many attempts have been made to assess the extent of the 
metallicity effect, with some analysing t he effect of Z on f i 
- a nd hence the PL z ero-point - such as iBono et al.l (|2008l ) 
and ISakai et al.l l|2004l ) , whilst others focus on the effect that 
Z has on the slope of the PL relation (e.g. lAlibert et al.l 
|l99p)). In this work we analyse the effect of Z on the de- 
rived distance modulus, which is quantified by the parameter 
7: 



7 = 



5{m - M)o 
SlogZ 



(4) 



where S{m — M)o is the difference between distance modu- 
lus corrected for the metallicity effect and without correc- 
tion, in the sense corrected - apparent, and SlogZ is the 



difference be tween a reference meta llicity and the region be- 
ing studied (i Kennicutt et al.] Il998l l . Recent measurements 
of 7 have produced negative values as large as 7 = —0.4 
mag dex~^, with the average being around 7 = ~0-25 mag 
dex ~^. Following the m ethodology of lKennicutt et al] (|l998t ) 
and lSakai et all (|2004h the value of 7 can be calculated using 
the two distance moduli, pinner and pouter, and the [0/H] 
gradient of M 3 3. We adopt the recent [O/H] gradient from 
iMagrini et all l|2007l ): 



d[0/H] 
dR 

d[0/H] _ 
dR 



= -0.19 ± 0.08 dex kpc"^ (R < 3 kpc) (5) 



-0.038 ± 0.015 dex kpc"^ (R > 3 kpc). (6) 



The deprojected radial distances of the centres of the two 
fields are 0.52 kpc and 5.11 kpc, leading to 12 -I- log (O/H) = 
8.852 and 8.296 for the inner and outer fields respectively. 
This gives a value of A[0/i/] = -0.556. To realise the fuU 
extent of the correction, we must initially treat the fields as 
if they had the same metallicity. This is done by using the 
same model PL relations for both regions. 

PL relations were fitted to the Cepheids in both fields in 
B,V,Ic an d the reddening-free Wesenheit magnitude Wvi as 



defined by lFouau6 et al.l l|2007l ) (F07), assumin g the Rv = 
3.23 reddening law described bv lCardeUi eTai] jl98a) 



W^, = V- 2.55(1/ - 7c). 



(7) 



The PL relations take the form M = a log P + b. The 
values for the slopes, a, were taken directly from the F07 
LMC relations, whilst the intercept, b, was derived from an 
iterative weighted least-squares fit for each relation. The val- 
ues of the Wvi zero-points are shown in Table |8l 
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Figure 6. Reddening-free Wesenlieit Wyi PL relations for the 
inner field (blue circles) and the outer field (red triangles). 



The zero-points of tlie PL relations for Wvi can be used 
to directly infer the distance modulus of M33, without the 
need to correct for extinction. The distance modulus can be 
found from the apparent Wvi magnitude of a Cepheid with a 
period of one day, such that jj, = b — Mo, where b is the zero- 
point of the PL relation and Mo is the absolute magnitude 
of a one day Cepheid, which is taken directly from the LMC 
PL relations from F07. 

The zero-points in Table |S] may at first glance appear 
discrepant. As they have been derived from the reddening- 
free Wvi PL relations we can safely rule out an erroneous 
extinction correction as the source of the difference. There 
are three other possible explanations for the difference we 
observe: 

(i) Blending affecting the magnitudes of Cepheids in the 
inner field. 

(ii) Differing period distributions of the two samples and 
a possible slope change at P = 10 days. 

(iii) Differing metallicity in the two fields. 



5.2.1 Blending 

The most obvious suggestion to explain discrepant values 
for the distance moduli is the fact the the inner field is ex- 
tremely crowded, and there may be Cepheids with unde- 
tected companions in the sample which will affect any PL 
relation zero-point that we try to fit. However, we do not 
believe that this is the cause of the observed deviation. 

Firstly, our procedure is based on allframe which uses 
a deep image, made by stacking many individual exposures, 
to do the photometry on each individual frame. When the 
photometry is done on this deep reference frame, the S/N 
for each object is much higher than on individual images, al- 
lowing us to better recognise, and classify as such, stars that 
may be b lended. Second l y, ou r comparison of the photom- 
etry with iMassev et al. I (|2006l ) shows excellent agreement. 



They use images with individual exposure times of 60 sec 
in BV and 150 sec in /; in other words they are up to 5 
times less deep than our data. This means that the level of 
crowding and blending will be less severe than in our case. 
Yet, our photometry is found to be in excellent agreement 
with theirs, for both fields, over the whole magnitude range 
covered by both surveys. If we were sensitive to blending 
levels that Massey et al. do not see, there should be a sys- 
tematic offset between the two catalogues. Ifowever, as we 
do not see a quantifiable offset in either Figure |4] or Table [7] 
we can safely rule out blending as the cause of the difference 
in distance modulus for the two fields. 

We can a lso rule out ble n ding as the cause on a the- 
oretical level. iFerrarese et al.l (|2000l ) addressed this issue, 
in the context of the HST Key-Project, with artificial star 
tests. They did this for 2 galaxies with fairly different levels 
of crowding. Their main conclusion was that this has little 
impact on distance moduli, provided the Cepheids used to 
derive the distance are carefully selected, as the ones in our 
final sample were. It is useful to note however, we found the 
same between our outer and inner fields using both the 
raw and selected Cepheid samples, although the PL rela- 
tion constructed from the selected sample has a much lower 
dispersion than that of the raw sample. This reduces the 
uncertainty on the fit, but also provides indirect evidence 
that blending is not an issue in our sample. 



5.2.2 Period Distribution and Slope Changes 

The second possibility for the differing values for the dis- 
tance modulus of M33 is the period range of the two sam- 
ples. The distribution of periods in the final, selected sample 
is shown in Figure [T] It is quite clear that the period dis- 
tribution of the two regions is not the same. This is for two 
reasons. Firstly, the short period Cepheids are fainter so we 
detect more of them in the outer field where the detection 
limit is lower. Secondly, long period Cepheids are rarer than 
their short period counterparts. Therefore, we expect the ra- 
tio of long period Cepheids to total number of stars detected 
to be quite low, hence they will be preferentially detected 
in the region with a higher number of stars; in this case the 
inner region. 

There is some observational evidence that the slope of 
the V and Ic PL relations change around 10 days (e.g. 
iNgeow et al] l|2009l )). However, there is still debate about 
whether this effect is real, and if the same non-linearity will 
appear in the Wvi relation. To rule out any possible slope 
change as the cause of the discrepant distance moduli the 
samples were both cut to logP < 1 and the PL relations 
refitted. Note that no attempt is made to determine if any 
non-linearity is present or to quantify the effect, just to elim- 
inate the possibility; such tests would require a much larger 
Cepheid sample than is used in this work. As a first test 
both the slope and zero-point of the fits were allowed to vary. 
The two slopes had the same value, as expected, but again 
the two zero-points were different. To bring the uncertainty 
on the zero-points down the slope was then fixed and the 
zero-points refitted. The resulting fit is shown in Figure H) 
Using the fixed slopes gives zero-points of 21.90 ± 0.02 and 
22.03 ± 0.04 for the inner and outer fields respectively. 

The fact that the zero-points are still discrepant at the 
3(T level when we use a fixed period range demonstrates that 
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Figure 5. PL relations for the inner field (left hand p anel) and the o uter fi eld (right hand panel), B,V and Ic go from left to right. The 
lines correspond to the LMC PL relations described in lFougue et al. I ||2007|) . The magnitudes have not been corrected for reddening 
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Figure 7. Period distribution of the two Cepheid samples used 
to fit the Wyi PL relations 
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Figure 8. W^i PL relation fit using logP < 1 criteria. The two 
slopes were found to be identical, however the zero-points still 
differ at the 3cr level. 



this effect, whatever the cause, is universal over the period 
range and is not due to the possible non-linearity of the PL 
relation. 



Table 8. Distance moduli from the W^i PL relations using the 
whole Cepheid sample and the period limited sample. 



5.2.3 Metallicity 

As we have successfully ruled out blending and the period 
distribution as causes of the discrepant distance modulus the 
only remaining explanation is the differing chemical abun- 
dances in the two environments. The zero-points and dis- 
tance moduli for each fit are shown in Table |S] 

To quantify the effect of metallicity we use 7, as de- 
fined in Equation U We can substitute S{m — M)o for 

— 0.16 ± 0.04 mag, using the full sample distance mod- 
uli, and SZ — — 0.556±0.101 dex. This leads to a metallicity 
correction of 7 = —0.29 ± 0.11 mag dex~^. Thi s is compat- 
ible w ith recent measurements, such as that of ISakai et al.l 
l|2004h . Our final distance measurement for M33 is calculated 
by assuming an LMC metalhcity of 8.34 jSakai et al.ll2004l ). 
Using 7 = — 0.29±0.11 we find that the metallicity corrected 
distance modulus to M33 is n-, — 24.53 ± 0.11, consistent 
with recent TRGB distance measurements. The individual 



Field bLMC fJ- 



All P 


Inner 


21.85 ±0.02 


24.37 ±0.02 


Outer 


22.01 ±0.03 


24.54 ±0.03 


log P < 1 


Inner 


21.90 ±0.02 


24.42 ±0.02 


Outer 


22.03 ±0.04 


24.55 ±0.04 



distance moduli derived for each Cepheid are shown in Fig- 
ure [1 

The uncertainty on our result is significant. By using 
a larger Cepheid sample the contribution from the distance 
moduli could be reduced. Also, the sample only covers the 
two extremes of the metallicity gradi ent. We will investigat e 
the effect further using the data from lHartman et al.ll|2006ll . 
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Figure 9. Measured distance moduli for each Cepheid as a func- 
tion of deprojected galactocentric distance and metallicity. The 
line is a weighted least squares fit to the individual data points. 



Their sample contains over 2000 Cepheids spread over the 
entire galaxy, and hence the entire [0/H] range. This will 
allow us to fully investigate the effect over a wide range of 
[O/H], rather than just fitting a straight line to two points, 
and the use of so many Cepheids will bring the random 
uncertainty down considerably. 

We note that the largest contribution to the error bud- 
get is the uncertainty on the metallicity gradient of M33. 
This could be r educed using a large sample of Cepheids. 
Recent work by iBeaulieu et al.l (|2006h has shown that the 
two periods of double mode Cepheids can be used to calcu- 
late their metallicity. By using a large sample of these stars 
across the whole galaxy, in combination with observations 
of fundamental mode Cepheids covering the same area it 
would be possible to calculate 7 with a much higher degree 
of accuracy. 



6 CONCLUSIONS 

We have produced a photometric catalogue in B,V and Ic 
for over 200,000 objects in two regions of M33, identifying 
167 fundamental mode Cepheids. Using the reddening free 
Wesenheit LMC PL relations we have established distances 
to the two regions as pinner = 24.37±0.02 mag and fiouter ~ 
24.53±0.03 mag. The difference in these two distance moduli 
is believed to come from the effect of metallicity on the zero- 
point of the PL relation, and hence its effect on Cepheid 
distances. We have quantified this effect by the parameter 

7 ~ ^^"1^ /^}^ ^^'^ fiiid a value of 7 = —0.29 ± 0.11 mag 

d[0/H\ 

dex~^, consistent with other recent estimates. For the range 
of metallicities discussed here it is appropriate to include 7 
in the PL relation. This leads to the numerical PL shown 
in Equation [8l assuming an LMC distance modulus of 18.4 



mag. The metallicity of the region in question is denoted by 

l0g[0/HUes. 



W^i = -3.320(±0.011) logP + 22.01(±0.03) 



-F0.29(±0.11)(log[O/J/] 1 



l0g[0/H]r 



(8) 



The analysis presented here assumes that only the zero- 
point of the PL relation changes with metallicity; the slope 
does not significantly change between the MW and LMC 
metallicities. We have repeated the analysis using the MW 
slope of the Wvi PL relation from F07 and we found the same 
value for 7. This is an important result as it means that 
wc can find reddening- and metallicity-corrected Cepheid 
distances for any population from two colour {V and Ic) 
photometry without the time consuming task of calculating 
internal Ay values for other galaxies. 

As noted above, the final result and particularly the er- 
ror are significantly affected by the value of metallicity one 
uses for each r egion . We have used t he recent results from 
iMagrini etlll (|20o3). ICrockett eTsU ifioOSl found a much 
shallower gradient however. Using their values of metallicity 
in our two fields, we find 7 = —2.899 ± 0.065 mag dex~^, 
incompatible with other determinations of this parameter. 
Such a large value of 7 would have dramatic effects on the 
distance scale. In particular, it would change the distance 
to the LMC by so much that the Cepheid distances would 
become incompatible with all other measurements. The dif- 
ference in distance moduli between the two fields is an argu- 
ment in favour of a rather large metallicity gradient in M33. 
However, the metallicity gra dient of M33 is s till su bject to 
strong debate and, to quote iKudritzki et al.l l|2008h . "...the 
expectation of homogeneous azimuthal metallicity in patchy 
star-forming s piral galaxies s e ems n aive...". Our extension 
of the work of IBeaulieu et all (|2006l ) based on double-mode 
Cepheids in M33 will shed light on this discrepancy by mea- 
suring the metallicity and its gradient using many more posi- 
tions. Hopefully this paper serves as an example of the need 
for better determinations of the metallicities of all galaxies 
where Cepheids have been observed. 

It is also important to note the correlation between th is 
work and the Hubble Key Project (jPreedman et al.ll200ll ): 
although the M33 Cepheids have shorter periods than the 
average in the HST project, the value of 7 that we find is 
consistent with their work. 

We will further investigate the effect of metallicity by 
applying th e methodology described above to the cata- 
logues from iHartman et al.l ([2006') . With the large number 
of Cepheids we will be able to sample the whole [O /H] range 
of the galaxy and the whole log P — Mag parameter space, 
finally calibrating the effect that has plagued the Cepheid 
PL relation and distance scale for many decades. 
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